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Abstract. 

We present the first L-band spectroscopic observations for a dozen stellar sources in the central 0.5 pc of the GC 

stellar cluster that are bright in the 2-4 /im wavelength domain. The L-band data were taken with ISAAC at 

the VLT UTl (Antu). With the aid of additional K-band spectroscopic data, we derive optical depth spectra of 

the sources after fitting their continuum emission with a single reddened blackbody continuum. We also derive 

^^ ' intrinsic source spectra by correcting the line of sight extinction via the optical depth spectrum of a late type star 

|/-\ I that is most likely not affected by local dust emission or extinction at the Galactic Center. The good agreement 

between the two approaches shows that the overall variation of the line-of-sight extinction across the central 0.5 pc 

V^ ' is AAk < 0.5 mag. The extinction corrected spectra of the hot He-stars are in good agreement with pure Rayleigh 

("^ I Jeans continuum spectra. The intrinsic spectra of all other sources are in agreement with continuum emission 

^T ■ and absorption features due to the dust in which they are embedded. We interprete both facts as evidence that 

^p ' a significant amount of the absorption takes place within the central parsec of the Galactic Center and is most 

(^ , likely associated with the individual sources there. We find absorption features at 3.0pm, SA/im, and 3A8fim 

O ,■ wavelength. Correlations between all three features show that they are very likely to arise in the ISM of the 

central 0.5 pc. Spacially highly variable hydrogen emission lines seen towards the individual sources give evidence 

of the complex density and temperature structure of the mini-spiral. In the cases of the sources with featureless 

K-band spectra like IRS 21 and IRS IW, they are consistent with massive hot stars embedded in the bow shock 

l— j I created by their motion through the dust and gas of the mini-spiral. It is likely that the bow shock scenario may 

be applicable to most of the dust embedded sources in the central stellar cluster. Spectroscopy of high MIR-excess 

sources 0.5" north of the IRS 13 complex is largely consistent with them being YSOs. However, a bow-shock 

J^> , nature of these sources cannot be excluded. The L-band spectrum at the location of SgrA*, is consistent with 

t . • that of a hot O-type star, such as S2, which was very close to Sgr A* at the time of our observations. 

Key words. Galaxy: center - galaxies: nuclei - infrared: stars - infrared: ISM extinction 

1. Introduction the tidal field of the black hole, intense stellar winds, and 

strong magnetic fields. 

At a distance of 8 kpc (Eisenhauer et al. 2003), the 

Near-infrared diffraction limited imaging over the past 10 Galactic Center is surrounded by a circumnuclear ring of 

years (Eckart & Genzel, 1996; Genzel et al., 1997; Ghez dense gas and dust showing clumpy extinction (Glisten et 

et al., 1998, 2000; Eckart et al., 2002; Schodel et al. 2002, al. 1987). Inside this ring, there is a central cavity of about 

2003; Ghez et al., 2003) has resulted in convincing evi- 3 pc diameter, that contains mainly ionzed or atomic gas. 

dence for a 3-4xlO^M0 black hole at the center of the The visual extinction estimates towards prominent sources 

Milky Way. This finding is supported by the discovery of within the central stellar cluster range between 20 and 

a variable X-ray and NIR source at the position of SgrA* 50 magnitudes with a median around 30 magnitudes (see 

(Baganoff et al., 2001; Genzel et al., 2003a). Most int- Rieke, Rieke, & Paul 1989, Chan et al. 1997, Scoville et 

rigingly, near-infrared imaging and spectroscopic obser- al. 2003). In addition Scoville et al. (2003) showed that 

vations have provided evidence for recent star formation the extinction is smoothly distributed across the central 

in the central parsec of the Milky Way, an environment 10 to 20 arcseconds with no indication of concentrations 

previously thought hostile for star formation because of of extinction on scales of about 1" to 2" . 
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The ~30 magnitudes of visual extinction along the line 
of sight toward the Galactic Center (GC) is mostly due to 
the diffuse interstellar medium (ISM) (Lebofsky 1979) and 
in part due to dense molecular gas (Gerakines et al. 1999; 
de Graauw et al. 1996; Lutz et al. 1996). The absorbing 
gas is cold (10 K) and the abundances of important molec- 
ular species are similar to the solar neighborhood (Moneti, 
Cernicharo, & Pardo 2001a, Chiar et al. 2000). In addi- 
tion Blum et al. (1996) and Clenet et al. (2001) concluded 
that the colors of individual dusty sources within the cen- 
tral stellar cluster contain a substantial contribution from 
intrinsic reddening. 

The entire central parsec of our Galaxy is powered by 
a cluster of young and massive stars (Blum et al. 1988, 
Krabbe et al. 1995, Genzel et al. 1996, Eckart et al. 1999, 
Clenet et al. 2001). Within that cluster the 7 most lumi- 
nous (L>105'^5 Lq ), moderately hot (T<10'*-5 K) blue 
supergiants contribute half of the ionizing luminosity of 
that region (Najarro et al. 1997, Krabbe et al. 1995, Blum 
et al. 1995). Such massive and hot stars were also found 
in dense clusters within the Galactic bulge, i.e. the Arches 
cluster (Cotera et al. 1992, see also Figer et al. 2002 and 
references therein) and the Quintuplet cluster (e.g. Figer 
et al. 1997). 

In addition to the massive blue supergiants, a popula- 
tion of dusty sources associated with bright dust emission 
can be found in the Galactic Center stellar cluster. After 
initial preceeding work by Becklin & Neugebauer (1968, 
1969) first individual mid-infrared sources in the central 
stellar cluster (among them IRS 1, 3 and others) were re- 
ported by Rieke & Low (1973) and Becklin & Neugebauer 
(1975). Later, IRS 1 was resolved into multiple compo- 
nents by Storey & Allen (1983), Rieke et al. (1989), Simon 
et al. (1990), and Herbst et al. (1993). Further high reso- 
lution imaging by ToUestrup et al. (1989) resolved IRS 6 
and IRS 12 into multiple components. 

In this paper, we discuss MIR sources that are located 
well within the central stellar cluster at projected dis- 
tances from Sgr A* of less than 0.5 pc (Fig. ^. Several 
sources like IRS 1, 3, and 21 are dominated by dust emis- 
sion and are strong at a wavelength of lOfim, whereas the 
supergiant IRS 7 is brightest at 2.2/xm. The nature of the 
dust enshrouded sources is still unknown. One among the 
best studied cases is IRS 21, which is strongly polarized 
(17% at 2fim ; Eckart et al. 1995, Ott et al. 1999, Krabbe 
et al. 1995). Initially, Gezari et al. (1985) suggested that 
IRS 21 is an externally heated, high-density dust clump. 
Given the MIR excess and the featureless NIR spectra 
several other classifications have been proposed, including 
an embedded early-type star and a protostar (Blum et al. 
1988, Krabbe et al. 1995, Genzel et al. 1996, Clenet et al. 
2001). Tanner et al. (2002) suggest that IRS 21 is an op- 
tically thick dust shell surrounding a mass-losing source, 
such as a dusty recently formed WC9 Wolf-Rayet star. 
Tanner et al. (2002, 2003) indicate that the extended dust 
emission of most of the central sources is consistent with 
bow shocks created by the motion of massive hot stars 
through the dust and gas of the mini-spiral. 



One way of investigating the nature of these bright 
NIR/MIR sources is imaging and spectroscopy in the 2 
to 4/ini wavelength range. In addition to Hydrogen and 
Helium recombination lines this wavelength domain is 
dominated by strong absorption features due to abundant 
molecules (NH3, CH3OH, H2O, CO, CO2 etc.), func- 
tional groups (like NH2, CH2), and ices. Here H2O ice 
enriched with molecular material is of special importance. 
Liquid, crystalline, amorphous water ice as well as trapped 
water ice in SiO condensate (Wada et al. 1991) give rise 
to a rich variety in shapes of a prominent feature with 
its deepest absorption at 2.94-3.00 ^m (e.g. Wada et al. 
1991). The variety in shapes of the water ice feature is 
dependent not only on temperature, but also on anneal- 
ing history and on the ice composition etc... (Hagen et 
al. 1983, Tielens & Hagen 1982, Tielns et al. 1983, Kitta 
& Kratschmer 1983, Hudgins et al. 1993, Maldoni et al. 
1998) 

The emission of dust and the absorption features of 
ices are important diagnostic tools for the investigation of 
the interstellar medium and circumstellar environments of 
individual sources. 

Infrared sources towards the Galactic Center show a 
wealth of ice absorption features (Butchart et al. 1986; 
Sandford et al. 1991) indicative for a broad range of or- 
ganic material mostly in the diffuse interstellar medium. 
Aliphatic hydrocarbons are characterized by their CH2 
(methylene) and CH3 (methyl) stretching modes around 
3.4/xm (Sandford et al. 1991; Pendleton et al. 1994). 
Aromatic hydrocarbons are detected via their CH and 
CC stretching modes at 3.28 and 6.2/Ltm. (Chiar et al. 
2000, Pendleton et al. 1994). An absorption feature at 
3.25/im has been found towards dense molecular clouds. 
It is attributed to aromatic hydrocarbon molecules at low 
temperatures (Sellgren et al. 1995; Brooke et al. 1999). 
Differences in the exact central wavelength and profile 
width of the absorption near 3.3/im are mostly attributed 
to differences in temperature and/or carrier of the absorb- 
ing molecules in these regions. 

In this paper we present 3 to 4/im imaging and spec- 
troscopy data combined with near- infrared 2.2/^ni spec- 
troscopy on the strongest mid-infrared sources in the cen- 
tral stellar cluster. In addition to the previously published 
L-band observations of IRS IW, IRS 3 and IRS 7 we 
provide the first L-band spectra of 9 other MIR sources: 
IRS 9, IRS 13, IRS 13N, IRS 21, IRS 29 and IRS 16 C, 
CC, NE and SW. These data on sources located in the 
central 0.5 pc of the GC enable us to study the proper- 
ties of the local interstellar medium and of circumstellar 
matter in this region. 

2. Observations and data reduction 

In order to investigate the nature of the dust embedded 
sources within the central 0.5 pc of the nuclear stellar 
cluster we used near- and mid-infrared imaging and spec- 
troscopy. In the following we describe the instrumentation 
and the data reduction that was employed. 
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2.1. MIR Observations and data reduction 

The mid-infrared imaging and spectroscopy was obtained 
using the ISAAC instrument at the ESO Very Large 
Telescope (VLT) unit telescope UTl (Antu), at the 
Paranal observatory in Chile. We have performed spec- 
troscopic and imaging observations of the Galactic Center 
during May 23-30, 2002, as part of a monitoring program 
of SgrA* (Eckart et al. 2003, Baganoffet al. 2003). AU the 
data reduction has been performed using routines from the 
IRAF and MIDAS software packages. 

For the imaging in the L-band, paired flat fielded im- 
ages at different chopper throws (18" ±2") and chopping 
position angles (0° to 180°) were subtracted from each 
other, resulting in frames containing a positive and a 
(shifted) negative image. The frames were then shifted 
to a common reference point that coincides with a pos- 
itive image of a source. Subsequently, frames belonging 
to the same batch, i.e. taken sequentially with identical 
or different chopper throws and/or chopping angles were 
combined by calculating a median. Since the images were 
moved to a common reference point this procedure elim- 
inates the negative "shadows" generated by the subtrac- 
tions. This procedure also effectively removes cosmic rays 
and bad pixel structures. Such a batch typically consists 
of up to 40 images, and each resulting combined image 
covers an integration time of approximately 35 minutes. 
For the present investigation we used the images with the 
best seeing and image quality with an angular resolution 
of about 0.4". The absolute L-band fiux calibration (see 
magnitudes listed in Table ^ was performed using the 
fluxes of several bright well isolated objects also measured 
by Clenet et al. (2001). 

The spectroscopic observations were performed with 
the long-wavelength (LWS3) and low resolution (LW) 
mode using the SL filter covering the wavelength range of 
2.7//m - 4.2/im and 4.4/im - 5.1/im, respectively. The use of 
a 0.6" slit width implied a spectral resolution [R = A/AA) 
of R=600 in that wavelength domain. The seeing at this 
time was ranging between 0.4 and 0.9 arcseconds. In 
order to compensate for the thermal background sepa- 
rate chopped observations were carried out using chopper 
throws of ^^18 arcseconds and random nodding within 2 
arcseconds along the slit. We adopted 4 different slit posi- 
tions all running through the SgrA* location (see Fig.QJ. 
The resulting images were divided by fiat-fields, corrected 
for cosmic rays, for sky lines and dispersion related dis- 
tortion. The wavelength calibration has been performed 
using a Xenon- Argon lamp. 

Two chopped frames (with shifted image positions) 
were then subtracted one from each other to provide a 
single frame containing two negative trace images and a 
positive one with twice the intensity than the negative 
images. After extraction of the individual source spectra, 
they were corrected for wavelength dependent sensitiv- 
ity, atmospheric transmission, and telluric lines using two 
standard stars HD 194636 (B4V) and HD 148703 (B2in- 
IV). 




Fig. 1. L-band image as obtained with ISAAC at the VLT 
UTl. The four slit positions chosen for the spectroscopic 
observations are also shown. They all run through the po- 
sition of SgrA*. The position of the CO-star used for the 
calibration along the line of sight is shown as well. 



The spectra were normalized to the given K- and L- 
band magnitudes listed in Tabled 



2.2. NIR Imaging spectroscopy and data reduction 

The near-infrared data in the K-band (2.2/im) was ob- 
tained via integral field spectroscopy using the imaging 
spectrograph 3D (Weitzel et al., 1996) combined with 
the tip/tilt corrector ROGUE (Thatte et al., 1995) at 
the ESO/MPG 2.2m located at La Silla. This instrument 
allows observations of a continuous 2-dimensional field 
(16 X 16 pixels) while providing spectral information for 
each spatial image element. These seeing-limited observa- 
tions were done with a pixel scale of 0.3" resulting in a 
field of view of 4.8" x 4.8" for individual pointings. Using 
a spectral resolution of R — 2000 each half of the K-band 
had to be covered separately. 

The observations were centered on the IRS 16 cluster. 
In March 1996, the central parsec was observed with a 
spectral coverage from 1.97 to 2.21/^m (lower half of the 
K-band). In total 17 different but overlapping regions were 
observed. The total sky area covered amounts to 16" in 
East- West and 10" in North-South direction. The upper 
half of the K-Band from 2.18 to 2.45/im was observed 
during a second observing run in April 1996. Here a total 
of 52 overlapping regions was observed. This resulted in 
an area of 25" in North-South direction and 20" in East- 
West. For further details of the observations see (Genzel 
et al. 1997, Ott et al. 2003). 
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In order to calibrate the wavelength scale, spectral 
lamp data (argon lamps in this case) were taken at the 
beginning or the end of each observing night. As a further 
calibration step, calibration sources with a known spec- 
trum were observed at a similar airmass as the Galactic 
Center. These standard stars were divided by a spectrum 
of the same stellar type (Kleinmann & Hall 1986) in or- 
der to remove stellar features resulting in an atmospheric 
transmission spectrum. The source data were then divided 
by this spectrum. 

3. Extinction correction 

The patchy extinction towards the GC stellar cluster (see 
e.g. Scoville et al. 2003) demand a careful cahbration of 
spectral data. Therefore, the correction of the extinction 
along the line of sight has been carried out using the 
method described in Sect. 13.11 and tested using another, 
independent method described in Sect. 13.21 A cross-check 
shows that both of our approaches agree very well (see the 
results in Sect. 13.2."^ and the conclusion fSect. l^^l) at the 
end of this section) and that the extraction of the optical 
depth spectra described in Sect. ^ can be done safely. 

3.1. Fitting Blackbody Emission to ttie 2 to 4fim 
Spectral Continua 

3.1.1. Method 

In order to estimate the extinction towards the individ- 
ual sources and to determine the approximate continuum 
shape of the observed spectra, we have performed simulta- 
neous fits of our K- and L-band data with single, reddened 
blackbody spectra. In this process we have considered the 
continuum emission, fitted locally around 2.2/ir7i, 3.75/xto, 
and 4.17/ir7i wavelength, as representing the intrinsic, red- 
dened continuum of the sources. In order to find the best 
fit, the temperature of the blackbody was allowed to vary 
within limits with a step size of 100 K. For the Helium 
stars these limits were 20,000 and 26,000 K (Najarro et al. 
1997). For the remaining objects the limits were 200 K and 
4,000 K. We allowed the K-band absorption Ak to vary 
over a range between 2.7 to 4.5 magnitudes with a step 
sizes of 0.05 mag assuming the extinction law A\ ^ A^^'^ 
stated by Martin & Whittet (1990). We extended the fit- 
ting range towards higher Ak values in order to allow a 
search for higher extinction towards the dusty sources. In 
general the fitting is not very sensitive with respect to the 
determination of Ak due to the strong H2O absorption 
feature that dominates almost the entire L-band. 

There are several potential problems that one has to 
consider when using this procedure (see also discussion in 
Chiar et al. 2002): 

We have to take into account that our K- and L-band 
data have not been taken at the same epochs and in some 
cases their relative calibration may be affected by vari- 
ability. Some sources in the central parsec of the Milky 
Way are variable on time scales of months to years (Ott 



et al. 1999, Blum et al. 1996, and Tamura et al. 1994, 
1996). In order to achieve a successful combined fit over 
both bands we therefore used different overall flux cali- 
brations of the spectra. For the K-band fluxes, we have 
considered data from 5 different references: Becklin et al. 
(1978), ToUestrup et al. (1989), Blum et al. (1996), Ott et 
al. (1999), and Clenet et al. (2001). The observations of 
Clenet et al. (2001) are the most recent of the four refer- 
ences and have been obtained with Adaptive Optics. For 
this reason, we have privileged the use of the values given 
in Clenet et al. (2001) for the fltting procedure. However, 
for some sources, a better flt was obtained when using 
a different flux calibration. This was the case for IRS 7, 
IRS 9, the stars in the IRS 16 complex as well as IRS 13. 
In the case of IRS 7 and IRS 9, these objects are vari- 
able as shown in Ott et al. (1999), Blum et al. (1996), 
and Tamura et al. (1994, 1996). The relative flux density 
calibration from Becklin et al. (1978) results in the best 
fits for our data since their observations in the K- and L- 
bands were performed in the same year. The IRS 13 com- 
plex consists of at least three bright members. El, E2, 
and E3 (Paumard et al. 2001, MaiUard et al. 2003). They 
contain the hottest and most luminous star in the entire 
region (Najarro et al. 1997). The IRS 13 complex therefore 
is very likely to be variable in flux density. In this case the 
fitting results were best using the K-band fluxes by Blum 
et al. (1996). 

The blackbody temperatures and K-band extinction re- 
sulting from the best fit of the spectra are listed in Tabled 

Moreover, most of the spectra are heavily affected by 
broad absorption features and it is difficult to determine a 
clear measure of the underlying continuum emission. Also, 
the assumption of a single temperature blackbody contin- 
uum can only be taken as a first approximation. On the 
other hand, reddened multi-temperature models quickly 
result in a larger number of not well determined parame- 
ters (e.g. temperature and relative flux density contribu- 
tion for each component). 

The resulting best fits of the individual spectra with a 
single blackbody are shown in Figs [3 to 0] 



3.1.2. Results 

We list the resulting parameters of the fits in Tableland 
the spectra with the corresponding blackbody curves in 
Figs. [3 13 and 01 The temperatures and the K-band ex- 
tinctions derived by the fitting procedure agree well with 
those found by other authors: The He stars of the IRS 16 
complex are best fitted with temperatures that are sim- 
ilar to those found by Najarro et al. (1997). IRS 7 and 
IRS 9 have typical temperatures of late-type supergiants 
or giant stars (Chiar et al. 2002, Ott et al 1999) and the 
remaining dusty sources are well fitted by blackbody con- 
tinua of typical temperatures for hot dust (Tanner et al. 
2002, Gezari et al. 1996, Genzel et al. 1997, Blum et al. 
1996 and others). 
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Fig. 2. Best fits of the K- and L- band spectra with single reddened blackbody continua. When no K-band spectrum 
was available, the flux density at 2.2fim represented by a diamond was considered in the fitting procedure. 
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Fig. 3. Fitting spectra with single reddened blackbody continua (cont.). 
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Fig. 4. Fitting spectra with single reddened blackbody continua (cont.). 



Source 


Blackbody 


Ak 


L-band 




Temperature 




magnitude 


IRS IW 


900 


3.30 


5.8 


IRS 3 


800 


4.20 


5.3 


IRS 7 


2100 


3.55 


5.0 


IRS 9 


1900 


3.10 


7.4 


IRS 21 


1200 


3.95 


6.8 


IRS 29 


1600 


3.70 


7.4 


IRS 16C 


22600 


3.45 


8.6 


IRS 16CC 


20200 


4.00 


9.0 


IRS 16NE 


22400 


3.25 


7.7 


IRS 16SW 


22500 


3.25 


8.4 


IRS 13 


1300 


3.00 


6.7 


IRS 13N 


1000 


3.90 


9.76 



Table 1. Results from fitting a single reddened blackbody to the K- and L-band spectra. In addition to the blackbody 
temperature and the K-band extinction obtained from the fitting procedure, we list the L-band magnitudes (±0.2 mag 
uncertainty) we derived from our images. 



The overall shapes we obtain by fitting single black- 
bodies are reasonable and compare favorably to those of 
Chiar et al. (2002) (especially in the case of the three 
sources in common IRS IW, IRS 3 and IRS 7). In ad- 
dition, they are supported by an independent procedure 
to calibrate the line of sight absorption described in the 
following section fSect. rO| . Also, the inclusion of the K- 
band spectra which was missing in the work by Chiar et 
al. (2002) allows to better judge the quality of the fit- 
ting procedure via comparison to the K-band continuum 



fluxes and spectral shapes, especially for the three sources 
in common (IRS IW, IRS 3 and IRS 7). 

Despite of the improvements in the spectral fitting, 
there are some mismatches we want to comment on: 
The fit of the IRS IW spectrum shown in Fig. [21 does not 
match the K-band spectrum of this source and shows a 
continuum level mismatch in the 3.8fim to A.Ofim region 
that has no physical significance. This is certainly due to 
a non-consistency between the flux calibrations in the K- 
and L- spectral bands. For this reason, we decided to set 
an upper and a lower limit to the fit by performing a fit 
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matching perfectly the K-band spectrum without caring 
about the adjustement at the red part of the spectrum 
on the one hand and another fit compensating the non 
physical mismatch at 3.8/im. 

An absorption at '^ 4.2/im appears in the spectra of 
IRS 9 and IRS 29. This suppression in flux is hardly an 
artefact of the data reduction, since we used the same 
procedures as for the other sources. We have no scientific 
explanation for that suppression but considered it as real. 
Also, both in shape and center wavelength, it does not 
seem to be consistent with the closest (at 4.27/xm) strong 
absorption line due to the stretching mode resonance of 
solid CO2 also observed towards the dust shells of some 
Young Stellar Objects (de Graauw et al. 1996). The pres- 
ence of this suppression does not affect results presented 
and conclusions drawn in this paper. 

Concerning the IRS 16 objects, all fits were not very 
satisfying towards the red part of the L-band spectrum 
and match perfectly well the K-band spectra. It is ob- 
vious that this is not due to the fitting procedure itself 
but to the shape of the L-band spectra of these sources 
which are very flat (except for the IRS 16 NE case where 
the shape matches better the overall shape of a Rayleigh- 
Jeans spectrum) . These spectra have been reduced in the 
same way as all the other spectra with the same flux cal- 
ibrator stars and thus the shape of the spectra is real. 
The origin of this behavior is not clear. It may be due 
to the fact that these sources are the least contaminated 
by local dust emission features and are likely most sus- 
ceptible to variations in the wavelength dependent line of 
sight extinction or of properties of the associated material. 
If for instance - compared to the intrinsic, local absorp- 
tions of the dusty sources - the line of sight absorption 
is dominated by amorphous H2O ice (Wada et al. 1991), 
the corresponding line center and strength of the red wing 
would be shifted towards the red. As a result the slope of 
the reddened continuum simultaneously fitting the K- and 
L-band continuum would be systematically too large. 

3.2. Determining the Wavelength-Dependent Line of 
Sight Extinction 

3.2.1. Method 

The calibration procedure described in this section is 
supported by the finding of Scoville et al. (2003) that the 
extinction shows a smooth distribution across the central 
10 to 20 arcseconds with no indication of concentrations 
of extinction on scales of about 1" to 2". Furthermore 
Blum et al. (1996) and Clenet et al. (2001) concluded that 
the colors of the sources within the central stellar cluster 
are due to both, intrinsic and foreground reddening. 

One of the slit settings we adopted for our obser- 
vations runs through a late-type star shown in Fig. ^ 
This star shows clear 2.3/im CO absorption band-heads 
in its K-band spectrum (see Fig. EJ, therefore, we will 
call it hereafter "CO-star" (e.g. Eckart et al. 1995). It is 



located at a projected distance of 12.6" (~ 0.5 pc) from 
the center, is well off the mini-spiral emission and does 
not show excess emission at wavelengths of 3/im or longer 
(see the L-band spectrum in Fig. Cj). 

Therefore we can safely assume that this star is largely 
free of local reddening and that its spectrum is mostly 
affected by the line of sight extinction along the 8 kpc 
towards the Galactic Center. 

We assume that the L-band spectrum of the CO-star 
can be represented by a 3600 K (M0-M3 spectral type) 
blackbody spectrum. If one matches the ~4.2/im flux den- 
sities of the measured L-band spectrum and the theoreti- 
cal continuum blackbody spectrum then the ratio between 
both spectra provides a measure of the relative wavelength 
dependent ice feature dominated L-band absorption due 
to the interstellar medium along the line of sight. The 
corresponding optical depth spectrum is shown in Fig. |S1 
If we correct this spectrum with the continuum optical 
depth value at 4.2/im of t=1.09 (derived via interpolation 
between the expected L- and M- band extinctions given by 
Rieke & Lebofsky 1985) the mean optical depth would be 
of the order of 1.64. This result is consistent with the value 
of 1.55 obtained by Rieke & Lebofsky (1985) for the L- 
band optical depth. We have divided the spectra of all re- 
maining objects by this extinction spectrum (of which the 
optical depth spectrum is shown in Fig.lHJ. The resulting 
spectra then represent probably more closely the source 
spectra as seen at the location of the Galactic Center - 
corrected for wavelength dependent absorption along the 
line of sight towards the Galactic Center. 

3.2.2. Results 

We show the spectra after correction for wavelength- 
dependent absorption together with our L-band image in 

Fig. El 

In Figs, ini to ^2 the spectra corrected for the measured 
wavelength-dependent line of sight extinction as deter- 
mined from the CO-star are shown in comparison to 
the continuum spectra that were derived by fitting the 
no77.-corrcctcd spectra with a single reddened blackbody 
rSect.rOTl. 

We find that they can be grouped into three classes rep- 
resented by the three figures. 

— I. As shown in Fig. ^2 the spectral shapes of the hot 
stars found in the IRS16 complex are fitted well by 
Rayleigh Jeans continuum spectra i.e. high tempera- 
ture blackbody spectra (listed in Table ^|. 
The slight mismatch between the blackbody spectra 
and the calibrated observed spectra in the red part 
of the L-band wavelength range corresponds to the 
one between the uncorrected spectra and the reddened 
blackbody continua shown in Fig.^ Its origin is, thus, 
probably the same as mentioned in the previous sec- 
tion. 
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Fig. 5. Our ISAAC L-band map combined with the spectra of the central sources corrected for the hne of sight 
extinction. 



II. The corrected spectra of IRS IW, IRS 3, IRS 13N 
and IRS 21 (Fig.^ are fiat in the red part with indi- 
cations of continuum absorption at wavelengths short- 
ward of 3.2/i7n. The spectra can be roughly explained 
by a T<1500 K blackbody continuum, with deviations 
at wavelengths < 3.5/xin due to the red wing of the 
3.0/xm ice absorption feature. 



III. Finally, in addition to strong H2O ice absorption 
the spectra of IRS 7, IRS 9, IRS 13, and IRS 29 (Fig. 
llOf) . show long- ward of 3.2/ir7i the characteristic shape 
of a T- 1500 K to 2000 K continuum. The absorp- 
tion at 4.2/iTO corresponds to that distinguished in the 
observed spectrum of Fig. |3 
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Fig. 6. The H- and K-band spectrum of the late-type CO- ^^g" '^- ^he L-band spectrum of the late-type CO-star 
star of which L-band spectrum is used for the line of sight ^^^d for the line of sight absorption correction, 
absorption correction. 



3.3. Conclusion on the extinction correction 

The comparison of the corrected spectra from extinction 
(using the CO-star) and the blackbody continuum spec- 
tra with temperatures equal to the ones obtained by the 
fitting procedure fSect. IXTl and Table 0) shows that the 
overall shape of the corrected spectra is consistent with the 
blackbody continuum shape of the corresponding balck- 
body temperature. The dusty sources are fitted well with 
blackbody temperatures of the order of ~ 1000 K and the 
global spectral shape of the hot Helium stars in the IRS16 
complex closely resembles pure Rayleigh Jeans spectra. 
This implies that the L-band optical depth spectrum of 
Fig. IHl is consistent with the extinction law of Martin & 
Whittet (1990) (see section ITT|l and with the known ab- 
sorption value 2.7 < Ak < 4 toward the Galactic Center 
(Rieke & Lebofsky 1985). Consequently, the fitting proce- 
dure made in Sect. 13. ll is reliable and one can use the red- 
dened continua obtained in that section in order to derive 
optical depth spectra. The extraction of the optical depth 
spectra is described in the next section. In the following, 
we do not make use of the corrected spectra obtained in 
Sect. 13.21 because small absorption features could still be 
present in the spectrum of the correction CO-star which 
may affect the optical depth measurments. 

In addition, the good agreement between the two cal- 
ibration procedures implies that the overall variation in 
extinction across the central 0.5pc of the Milky Way can- 
not be much larger than AA/^ =±0.5 magnitudes. This is 
consistent with the results by Scoville et al. (2003) who 
derived extinction estimates from the Pa/6cm radio con- 
tinuum and the Pa/H92a line emission over this central 
region. These estimates result in a distribution which is 
smooth on the ~1" scale. This supports the assumption 
that the excess extinction seen towards some of the sources 
must in fact be associated with the individual objects 
rather than with the diffuse ISM (see Blum et al., 1996; 
Clenet et al., 2001). 




^{|l m) 



Fig. 8. Optical depth spectrum of the line of sight extinc- 
tion towards the Galactic Center. 

4. Optical Depth Spectra 

4.1. Extraction of Optical Depth Spectra 

As it has been pointed out in the previous section, the 
optical depth spectra of the Galactic Center sources can 
be obtained using the reddened blackbody continua of 
the fitting procedure fSect. lTTI and Figs.|21to0J. 

With the results of the blackbody fitting (Table ^1, we 
have derived the optical depth spectra of the sources us- 
ing the equation Fobs = Fintre^~ where r is the optical 
depth, and Fobs and Fintr are the observed and intrinsic 
fluxes, respectively. All the spectra are shown in Fig. IT^ 
where they have been normalised to unity at 3.0/ir7i. Mean 
optical depth spectra for the dust enshrouded sources and 
the non-obscured sources are shown in Fig. 1131 This fig- 
ure shows very clearly three absorption features at S.Ofim, 
3.4/X7TI and 3.48/im. From Fig.^J it is evident that the ab- 
sorption depths of the 3.4/xm, 3.48fim and 3.0/im features 
vary from source to source. 

Mean spectra of the hydrocarbon feature at S.Afim and 
the ice feature at 3.0fim have been constructed in the same 
iterative way as in Chiar et al. (2002). The optical depth 
spectrum of IRS 29 shows the smallest 3.4/iTO absorption 
in comparison to the 3.0/ir7i feature. A first mean spec- 
trum of the S.Afim feature was obtained by subtracting 
the optical depth spectrum of IRS 29 from all other spec- 
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Fig. 9. L-band spectra corrected for line of sight absorption using the optical depth spectrum shown in Fig. |H1 
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Fig. 10. L-band spectra corrected for line of sight absorption (cont.). 
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Fig. 11. L-band spectra corrected for line of sight absorption (cont.). 



tra. Then the average spectrum of the residual spectra 
was derived. The mean water profile spectrum is then ob- 
tained by subtracting the average spectrum of the SAfim 
hydrocarbon feature multiplied by a free scaling param- 
eter. This parameter is adjusted such that a featureless 
spectrum is obtained in the SAfim region after subtraction 
from the mean source spectrum. Finally, the spectrum of 
the SAfim hydrocarbon feature was deduced by subtract- 
ing the mean S.Ofim spectrum from all the spectra and 
averaging the residuals. The resulting spectra are shown 
in Figs. [TH and [T31 

The most reliable results were obtained for the dust 
enshrouded sources. As a consequence of the successful 
continuum fits described in Sect. l3.il the absorption spec- 
tra level out at r^^O.O for wavelengths well above and well 
below the corresponding absorption features. For the hot 
He-stars in the IRS 16 cluster the situation is different. 
Here the continuum fits were less successful (see discussion 
in Sect. 13. l|) and the derived profile of the absorption fea- 
tures are less reliable towards the long wavelength end of 
the spectra. This is less severe for the 3.0fim feature since 
in this case only the continuum long-ward of '^S.S/zm is 
affected. 



4.2. Characteristics of Features in the Optical Depth 
Spectra 

An examination of the optical depth spectra reveals in- 
formation on the origin and physical conditions of the ab- 



sorbing material in the central parsec of the GC stellar 
cluster. We distinguish the following features: 

— The 3.0/i?n feature: The dusty sources (which are lo- 
cated in the minispiral area like IRS IW, 3, 13, 13N 
and 29) and hot He-stars (i.e. showing He emission 
lines in their K-band spectra which is the case for the 
IRS 16 sources, e.g. Najarro et al. 1997) allow to de- 
rive both the linewidth and line minimum of the 3.0/xm 
absorption feature. Both quantities are not affected by 
the slight continuum mismatch at wavelengths above 
3.5/im in case of the He-stars. The S.Ofxm absorption 
feature is associated with interstellar H2O ice in molec- 
ular clouds (Gillett & Forrest 1973). Laboratory ex- 
periments suggest that the shape of the observed ice 
absorption (which is deepest at about 2.95/im in IRS 7 
for example) toward the Galactic Center corresponds 
to trapped water ice in SiO condensate (Wada et al. 
1991). 

Our mean absorption spectrum (Fig. 114(1 shows, how- 
ever, that in our observations the deepest absorption 
occurs ^ 0.05/im longward of 3.0/im. This behavior 
is consistent with the presence of amorphous H2O ice 
(Wadaet al. 1991). 

Our mean spectrum of the 3.0/im ice feature appears 
to be narrower than the one obtained by Chiar et al. 
(2002). As modeled in their paper, the 3.0/im ice fea- 
ture derived by the authors is best fitted by a temper- 
ature of lOK and a maximum mantle thickness of 0.75 
to 0.85/Lim assuming the simple model of core-mantle 
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Fig. 12. Optical depths of the Galactic Center sources 
normalized to unity at 3.0/iTO. 



grains of Bohren & Huffman (1983) and the hypothesis 
of a variable mantle thickness of Smith et al. (1993). 
The S.^^m and 3.48fim features are known to arise 
in the diffuse interstellar medium. They are both part 
of a single line complex (Duley & Williams 1983) and 
have already been seen towards many Galactic Center 
sources (Sandford et al. 1991, Pendleton et al. 1994, 
Chiar et al. 2002, Mennella et al. 2003). They are most 
likely due to aliphatic hydrocarbons which are charac- 
terized by their CH2 (methylene) and CH3 (methyl) 
stretching modes (Sandford et al. 1991, Sellgren et 
al. 1995, Brooke, Sellgren & Geballe 1999, Grishko & 
Duley 2002, Butchart et al. 1986, Duley & Wilhams 
1984). 

Actually, the 3.4/im feature is composed of two bands, 
one at 3.38/im and the other at 3.42/im. These are due 
to the asymmetric stretching vibrations of CH^ and 
CH2 groups respectively (Sandford et al. 1991). 
The 3.48^m seen in all spectra (see also the mean opti- 
cal depth spectrum in Fig. lT^ is due to the perturbed 
symmetric C-H stretching vibrations of the same CH^ 
and CH2 groups (Sandford et al. 1991). This fea- 
ture was already reported by Allen & Wickramasinghe 
(1981), Jones et al. (1983) and Butchart et al. (1986). 
The 3.4^m absorption feature for the dusty sources 
and the He-stars are shown in Fig. El For both types of 
sources the short wavelength wing is very similar. The 
long wavelength wing of the line feature derived from 
the He-stars is affected by the continuum mismatch 
discussed in Sect. 13.11 



4.3. Measurements of the strengths of the Absorption 
Features 

The optical depths of the 4 absorption features at 3.0/iTO, 
3.4/i7TT, and 3.48/iTO are listed in Tabled The values are 
obtained by averaging the optical depth spectra over the 
wavelength intervals ]2.95/i7Tj,3.05/^7TT.[, ]3.35/iTO,3.45/i7TT.[ 
(including thus the two bands of the feature) and 




3.0 (im H2O 



3.5 
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Fig. 13. Mean L-band normalized spectra. The black and 
grey (or red in a colour version) spectra are obtained 
by averaging the dusty sources and the stellar (He-stars, 
IRS 7 and IRS 9) optical depth spectra, respectively. 
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Fig. 14. Mean isolated and normalized ice spectra of the 
dusty sources (in black) and He-stars, IRS 7 and IRS 9 
sources (in grey or red in a colour version). 



]3.43/iTO,3.53/im[, respectively. The 3.4/iTO and 3.48/1777 op- 
tical depths are measured in the spectra from which the 
3fim ice absorption feature was subtracted (see Sect. 14. 1|) . 
The error bars listed in the table take only into account 
the measuring uncertainties, except for the 3.0fim feature 
of the IRS IW case where the error bar is taken such that 
the optical depth value includes its upper and lower limits 
obtained by fitting upper and lower baselines as described 
in Sect. 1X0 

To estimate the lower limit of the optical depth values 
of the 3.4/ir77 absorption feature, we measured the value 
of the strength of this feature using a linear continuum 
between 3.25fj,m and 3.65/i77i as baseline. The obtained 
values are listed in Table [5] 

Table |21 shows that the optical depth values span a 
large interval suggesting that part of the absorption fea- 
tures arise probably from the local medium and may be 
associated with the individual sources. 

Concerning the optical depth values of the S.Ofim fea- 
ture, they agree well with the values obtained by Chiar 
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Object 


Ta.O/jm 


T3.4nm 


T3.4Aim lower limit 


T3.48Mm 


IRS13 


0.44 ± 0.03 


0.26 ± 0.02 


0.16 ±0.02 


0.19 ±0.02 


IRS13N 


0.35 ± 0.04 


0.33 ± 0.04 


0.21 ±0.03 


0.25 ± 0.03 


IRS16C 


0.67 ±0.04 


0.49 ± 0.05 


0.14 ±0.02 


0.45 ± 0.03 


IRS16CC 


0.425 ±0.05 


0.36 ± 0.03 


0.16 ±0.03 


0.31 ±0.02 


IRS16NE 


0.38 ±0.05 


0.27 ±0.02 


0.14 ±0.02 


0.23 ± 0.02 


IRS16SW 


0.57 ±0.03 


0.36 ± 0.05 


0.12 ±0.02 


0.29 ± 0.03 


IRSIW 


0.51 ±0.28 


0.32 ± 0.04 


0.16 ±0.02 


0.27 ±0.02 


IRS21 


0.58 ±0.03 


0.23 ± 0.04 


0.15 ±0.03 


0.16 ±0.02 


IRS29 


1.12 ±0.09 


0.37 ±0.03 


0.27 ±0.03 


0.29 ± 0.02 


IRS3 


0.30 ± 0.02 


0.23 ± 0.03 


0.20 ± 0.02 


0.15 ±0.01 


IRS7 


0.68 ±0.03 


0.41 ±0.06 


0.20 ±0.05 


0.33 ± 0.04 


IRS9 


0.67 ±0.06 


0.27 ±0.07 


0.16 ±0.04 


0.19 ±0.03 



Table 2. The optical depths of the 4 absorption features at S-Ofim, 3.3fim, SAfim and 3.48/im,. 
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Fig. 15. Mean isolated and normalized 3.4/iTO feature of 
the dusty sources (in black) and He-stars, IRS 7 and IRS 9 
sources (in grey or red in a colour version). The location 
of the 3.3/X7TI feature is indicated. 




3.4 

X(p,m) 

Fig. 16. Mean isolated and normalized 3.4^m feature of 
the dusty sources. The location of the 3.3fim feature is 
indicated. 



et al. (2002) for the two sources in common IRS IW and 
IRS 7. The 3.4/ir7i feature of these two sources exhibits 
higher values than those of Chiar et al. (2002). This is 
probably the result of the derivation of the isolated ice 
absorption feature described in Sect. 14.11 which is sensi- 



tive to the choice of the subtracted spectrum (here the 
spectrum of IRS 29). In particular, in Chiar et al. (2002), 
the authors find a negative optical depth value at 3.3^?7Z 
for IRS 7 due to the shape of the ice feature which is 
narrower than the mean ice absorption feature. Moreover, 
their optical depth value of the 3.4/zm feature of IRS 7 
is even smaller than the lower limit obtained in our spec- 
trum. 

Concerning the IRS 3 source results, no comparison can 
be done with Chiar et al. (2002) as the fit of their spec- 
trum was not satisfactory and did not match at all the 
K-band point flux. This may be due to the problematic 
relative flux calibration of their L-band spectrum which 
shows a different slope than ours. 

On the other hand, the lower limits of the optical 
depth values of the hydrocarbon feature for the IRS 3 
and IRS 7 sources in common with Sandford et al. (1991, 
1995) and Pendelton et al. (1994) agree well with those 
obtained by these authors. This is very satifactory as 
the values provided in these papers were derived using a 
linear fit similar to the one used here to derive the lower 
limits. 



4.4. Correlations Between the Strengths of the 
Absorption Features 

The plots of the optical depths of the absorption fea- 
tures at 3.4/X771 versus 3.0/im and i.4fim versus 3.48^to 
are shown in Figs. 1191 and 1181 respectivelv. In these figures 
we also plot the data by Chiar et al. (2001) that were 
mostly obtained on sources in the outer parts of the cen- 
tral stellar cluster. 

We find a good correlation (correlation coefficient of 
^ 0.98) between the 3.4/im and 3.48/im absorption fea- 
tures (Fig. I18|l and a trend of correlation (correlation co- 
efficient of ^ 0.43) between the 3.4/i77T, and S.Ofim features 
(Fig. I19|l . In Figs. 1181 and 1191 the best fits of linear regres- 
sion are drawn as well; the goodness of fit probabillities 
are respectively 0.9996 and 8 10~*. The first correlation 
is forseeable since both features (at 3.4^m and 3.48^to) 
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Fig. 17. L-band intrinsic spectrum (i.e. corrected from ex- 
tinction along the line of sight as explained in Sect. l3.2|l at 
the position of SgrA* taken through a 0.6" wide slit using 
ISAAC at the VLT UTl. The continuum of a blackbody 
of 22,200 K is also shown. The locations of hydrogen lines 
are marked. 



arise from the same functional groups as explained in Sect. 

113 

The possible correlation between the optical depths of 
the 3Afim and 3.0fim features shown in Fig. 1191 suggests 
that the ISM in the central region is a mixture of diffuse 
and dense material. On the other hand, if this correlation 
is real, then in this plot the discrepancy between the po- 
sitions of IRS 7 derived by us and by Chiar et al. (2002) 
cannot only be due to the variability of the source. In the 
case of variability, both positions should follow the overall 
trend of the correlation. Actually, the discrepancy between 
the positions of the 3 sources (IRS IW, IRS 3 and IRS 7) 
in common with Chiar et al. (2002) is due to the same 
reasons evoked previously in Sect. 14.31 

The fact that the 3.4/iTO to 3.0/ir7i absorption feature 
"correlation" shows a significant offset from the origin of 
the coordinate system used for display may indicate that 
there could be a significant SAfxm line absorption even if 
the water ice feature were not present. Therefore we may 
conclude that a certain amount of the 3.4/im line absorp- 
tion is due to the diffuse ISM on the line of sight to the 
Galactic Center. Due to the low density of the mini-spiral 
gas a substantial amount may, in fact, be closely linked to 
the individual sources that we are studying here in detail. 
Therefore, the remaining portion of the SAfim absorption 
may occur in the central parsec. This is also supported by 
the finding of Sandford et al. (1995) who derived anoma- 
lously high aliphatic CH absorption per visual extinction 
relative to all the other lines of sight for which data was 
then available. 



4.5. Hydrogen emission lines 

The obtained L-band spectra also allow us to investi- 
gate the HI line emission towards selected lines of sight. 
The corresponding line strengths, widths and equivalent 
widths for Bra, Pf7, and Pi6 are listed in Tables 01 
O and El At the achieved spectral resolution we can- 
not distinguish between the emission from the individual 
objects and the mini-spiral. The latter seems to domi- 
nate the emission since sources on or close to the mini- 
spiral (IRS IW, 16NE, 16SW, 21, 13, 13N), also show the 
strongest line emission. Especially for IRS 21 the K-band 
spectrum by Ott et al. (1999) suggests that probably all 
the hydrogen line emission is due to the mini-spiral. 
Using the Br^ and Pa line strengths published by Najarro 
et al. (1997) we can calculate the diagnostic line ratios 
Br^/Bra and Pa/Bra for three of the IRS 16 sources. 
The results are shown in TableOl In order to obtain dered- 
dened line strengths, we have considered an absorption in 
the L-band wavelength range of Al = 1.68 consistent with 
our optical depth spectrum of Fig. |S1 as well as the value 
given by Rieke & Lebofsky (1985). 

Taking the dereddened ratios at face value and omit- 
ting possible small contributions from the He line emis- 
sion close to Pa and Br7 line we can derive first esti- 
mates about the temperature and density of the emit- 
ting gas. Assuming case B recombination, for both sources 
IRS 16NE and IRS 16C the line ratios are consistent with 
emission from gas at low temperatures of < 5000K and low 
densities of <10^cm~'^. For IRS 16SW the ratios are more 
consistent with emission from a gas at high temperature 
~ 20, OOOK and high density ~ lO^'cm^'^. While it cannot 
be excluded that the emission towards IRS 16SW contains 
a significant contribution from the hot stellar atmospheres 
of the He-star IRS 16SW and the other neighboring hot 
stars in the IRS 16 complex, this result reflects the com- 
plex density and temperature structure within the mini- 
spiral that consists of a thermal plasma of Ue^lO'^cm"^ at 
T~10^ K (e.g. Brown, Johnston, & Lo 1981) with denser 
entities of n//>10^cm~^ and a few 100 K (e.g. Jackson et 
al. 1993). 

4.6. L-band spectrum at the position of SgrA* 

Our ISAAC L-band imaging and spectroscopy was part 
of a flux density monitoring program of SgrA* (Baganoff 
et al. 2001, 2003 and Eckart et al. 2004) which has been 
carried out simultaneously to CHANDRA observations. 
Therefore all slit settings were defined such that the po- 
sition towards SgrA* always fell into the slit. A mean 
extinction corrected spectrum at that position is shown 
in Fig. 1171 Here the weak hydrogen recombination lines 
are most likely due to emission from the mini-spiral. The 
Rayleigh Jeans shape of the continuum spectrum com- 
pares well to that of other hot stars, like the He stars in the 
central cluster. In Fig. El we show the L-band spectrum 
of SgrA* compared to a 22,200K blackbody spectrum. In 
fact, high resolution adaptive optics L-band images were 
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Object 


(P. + He)/Br^ 


{Br^ + He)/Br^ 


IRS16NE 

IRS16C 

IRS16SW 


3.2 ±0.6 
3.4 ±0.6 
4.6 ±0.8 


0.28 ±0.03 
0.26 ±0.03 
0.43 ±0.04 



Table 3. Diagnostic line ratios for the three IRS 16 sources in common with Najarro et al. (1997). 
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Fig. 18. Optical depth of the hydrocarbon absorption at 
the 3.4/zm versus the 3.48^m absorption feature. 
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Fig. 19. Optical depth of the hydrocarbon absorption at 
3.4/iTO versus the depth in the 3.0/x?n absorption fea- 
ture for all observed Galactic Center sources (the IRS 16 
sources are represented by stars and the remaining sources 
by triangles). The values obtained by Chiar et al. (2002) 
are also shown (in circles) but are not considered in the 
linear regression fit. 



taken in August 2002 with NAOS/CONICA on the VLT 
UT4 during the science verification phase (Clenet et al. 
2003, Genzel et al. 2003, Eckart ct al. 2004). These images 
show that at the position of SgrA* the L-band flux den- 
sity is dominated by the fast moving star S2 (Schodel et 
al. 2003). Our spectrmn is consistent with the fact that 
based on its K-band luminosity spectrum S2 is likely to 
be a ^15-20 solar mass late O, early B main-sequence star 
of age less than 20 Myr (Gezari et al., 2002; Eckart et al., 
1999; Figer et al., 2001; Ghez et al., 2003). 



5. IRS 13N 

Well within the central stellar cluster of the Milky Way, 
about 0.5" north of the IRS 13 complex,trend in Eckart et 
al. (2004) have found a small cluster (<0.13 hght years) 
of compact sources with a strong infrared excess and ob- 
tained a first spectrum of it. The integrated spectrum of 
that area clearly shows that this excess is due to the con- 
tribution of warm T=500 - 1000 K dust. The nature of 
these newly found sources is currently unclear. Eckart et 
al. (2004) discuss three possible explanations. The L-band 
excess sources north of the IRS 13 complex could be heav- 
ily extincted luminous stars. They could also be hot stars 
that heat and interact with the more ambient environment 
of the local mini-spiral. Finally these newly found objects 
may even be young stars and their L-band excess may be 
due to the flux density contribution of luminous accretion 
disks. While a combination of the first two possibilities 
currently cannot be entirely excluded, the luminosities and 
colors of these sources are consistent with those of young 
stellar objects. As described in Sect. 13. ll a reddened black- 
body continuum was also fitted to the IRS 13N spectrum 
considering the K- and L-band fluxes given by Eckart et 
al. (2004) (see Figs. El and EHl and Table [TJ. The result- 
ing temperature and K-band extinction agree well with 
the above explanations. The extinction of about 3.9 mag 
in the K-band is by about Imag higher than the overall 
line of sight extinction toward the Galactic Center and a 
~1000iir dust temperature emission agrees well with the 
highly dust embedded YSO's as studied by e.g. Ishii et 
al. (1998). The ice absorption feature often observed to- 
wards YSO's (see Ishii et al 1998) is not detectable in the 
extinction corrected spectrum shown in Fig. (SJ this may 
in part be due to the small absorption in the spectrum of 
the late-type CO-star that reverberates on the spectrum 
of the extinction along the line of sight (Fig. |Hl . It is also 
possible that the circumstellar material is affected by the 
intense radiation fleld at the Galactic Center (Lutz et al. 
1996). In this case the ices in the disks of young stars 
may become either tenuous and/or destroyed by hard-UV 
photons. 



6. Summary and Conclusions 

Combined NIR and MIR spectroscopic observations of 
sources in the central 0.5 pc parsec of the Milky Way al- 
lowed us to obtain a detailed picture of the absorption 
features visible in the spectra towards that region. Our 
investigation on the central sources (radius <2.5") com- 
plements the study by Chiar et al. (2002) that mostly 
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Object 


Line-Continuum of HI (4.051/im) 


Integral of HI (4.05 l^m) 
10-1* Wcm-^ 


EW of HI (4.051/im) 

A 


Line width 

A 


IRS 1 W 


55.9 ±0.40 


0.5869 ±0.1216 


-41.48 ±1.86 


210 ±42 


IRS 3 
IRS 7 


5.80 ±0.50 


0.0855 ± 0.0378 


-6.02 ± 2.83 


295 ± 105 


IRS 9 


19.60 ±0.40 


0.1137 ±0.0121 


-63.56 ±1.64 


116 ± 10 


IRS 21 


13.00 ± 0.20 


0.1566 ±0.0362 


-48.03 ±1.90 


241 ± 52 


IRS 29 


1.70 ±0.20 


0.0196 ±0.0041 


-10.26 ±1.75 


231 ±21 


IRS 13 


47.50 ± 0.40 


0.8479 ± 0.0570 


-87.59 ±2.63 


357 ±21 


IRS 13 N 


35.60 ± 0.30 


0.2243 ± 0.0393 


-63.58 ±1.18 


126 ±21 


IRS 16 CC 


1.69 ±0.03 


0.0275 ± 0.0066 


-58.20 ±4.39 


326 ± 73 


IRS 16 C 


3.86 ± 0.06 


0.0588 ± 0.0028 


-117.43 ±3.28 


305 ± 10 


IRS 16 NE 


10.19 ±0.08 


0.1126 ±0.0065 


-112.82 ±2.94 


221 ±11 


IRS 16 SW 


6.29 ± 0.07 


0.1025 ±0.0109 


-120.66 ±3.61 


326 ± 31 



Table 4. Continuum and emission line strengths as well as linewidths and equivalent widths of the Bra line towards 
the observed sources. 



Object 


Line-Continuum of HI (3.739/im) 


Integral of HI (3.739/im) 
10-1** Wcm.-^ 


EW of HI (3.739/im) 

A 


Line width 

A 


IRS 1 W 


13.50 ±0.60 


0.1276 ±0.0198 


-12.86 ±1.36 


189 ±21 


IRS 3 


3.30 ± 0.40 


0.0346 ± 0.0077 


-4.79 ± 1.06 


210 ±21 


IRS 7 


18.00 ±1.00 


0.1224 ±0.0158 


-4.41 ± 0.52 


136 ±10 


IRS 9 


3.80 ± 0.70 


0.0380 ± 0.0091 


-12.63 ±1.67 


200 ±11 


IRS 21 


2.10 ±0.10 


0.0166 ±0.0041 


-9.29 ±0.69 


158 ± 32 


IRS 29 


1.10 ±0.10 


0.0072 ± 0.0009 


-3.73 ±0.31 


131 ±5 


IRS 13 


8.90 ± 0.25 


0.0703 ± 0.0064 


-15.29 ±0.86 


158 ±10 


IRS 13 N 


6.30 ± 0.30 


0.0627 ±0.0127 


-15.38 ±1.64 


199 ±31 


IRS 16 CC 


0.37 ±0.03 


0.0039 ± 0.0007 


-9.96 ±1.14 


210 ±21 


IRS 16 C 


0.86 ± 0.05 


0.0086 ± 0.0009 


-20.03 ±1.84 


200 ± 10 


IRS 16 NE 


2.16 ±0.08 


0.0239 ± 0.0042 


-21.74 ±2.50 


221 ±31 


IRS 16 SW 


1.07 ±0.05 


0.0096 ±0.0010 


-19.72 ± 1.53 


179 ± 10 



Table 5. Continuum and emission line strengths as well as linewidths and equivalent widths of the Pf7 line towards 
the observed sources. 



includes sources with larger separations from the center. 
We find some evidence that the diffuse ISM in the 0.5 pc 
has properties that are slightly distinct from the ISM at 
larger distances from the center. 

The 3.0/im ice profile usually observed toward the 
Galactic Center peaks at 2.96/im, short-ward of the ice 
feature in local molecular clouds (McFadzean et al. 1989; 
Tielens et al. 1996; Chiar et al. 2000). It is likely due to 
cold (15 K) water ices with an enhanced NH3 abundance 
(Chiar et al. 2000). Especially processed ices contribute a 
substantial portion of the refractory grain materials that 
persist when the molecular cloud is dispersed by star for- 
mation, and these products may yield the extinction char- 
acteristic of the diffuse interstellar medium. 

However, the profile of the S.O/im ice feature obtained 
in our work peaks long- ward of S.Ofim and therefore is 
likely to be associated with the presence of amorphous 
H2O ice towards these sources. 

Simultaneous fits of our K- and L-band spectra with 
single reddened blackbody continua allowed us to estimate 
the extinction towards the individual Galactic Center 
sources and to determine the approximate continuum 



shape of the observed spectra. The derived K-band ex- 
tinctions and blackbody temperatures are consistent with 
values found for the sources in the central stellar cluster 
(Chiar et al. 2002, Tanner et al. 2002, 2003). 

Using the spectrum of a late-type star assumed to be 
free from local extinction, we were able to derive the spec- 
trum of the L-band extinction along the line of sight to- 
ward the Galactic Center. This spectrum has been used 
in order to derive the extinction corrected spectra of the 
most luminous sources. The extinction corrected spectra 
are consistent with the blackbody temperatures derived 
from the previous fitting procedure. The Rayleigh Jeans 
continuum spectra obtained towards all the hot stars in- 
dicate that the distribution of extinction in the central 
half parsec is fairly flat and varies in the K-band only 
by AAj^:~ ±0.5 mag. This is consistent with the results 
by Scoville et al. (2003). Therefore the excess extinction 
that we determined after correction for the line of sight 
extinction towards some of the sources should be associ- 
ated with the individual objects or very clumpy features 
of the ISM. This is in support of the findings by Blum et 
al. (1996) and Clenet et al. (2001) and is also consistent 
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Object 


Line-Continuum of HI (3.296/im) 


Integral of HI (3.296^m) 
10-18 Wcm-^ 


EW of HI (3.296Mm) 

A 


Line width 

A 


IRS 1 W 
IRS 3 


7.80 ± 0.20 


0.0573 ± 0.0097 


-9.26 ±0.81 


147 ±21 


IRS 7 
IRS 9 
IRS 21 


12.00 ±2.00 


0.0450 ±0.0135 


-3.50 ± 0.94 


75 ±10 


IRS 29 


— 


— 


— 


- 


IRS 13 


4.80 ± 0.20 


0.0226 ± 0.0033 


-7.65 ± 0.68 


94 ±10 


IRS 13 N 


3.30 ± 0.30 


0.0208 ± 0.0054 


-8.29 ± 0.99 


126 ±21 


IRS 16 CC 


- 


- 


- 


- 


IRS 16 C 


0.57 ±0.12 


0.0027 ±0.0008 


-7.57 ± 1.19 


94 ±10 


IRS 16 NE 


1.15 ±0.16 


0.0084 ± 0.0024 


-7.62 ±1.31 


147 ±21 


IRS 16 SW 


0.38 ± 0.09 


0.0026 ± 0.0008 


-4.30 ±1.18 


137 ±10 



Table 6. Continuum and emission line strengths as well as lincwidtlis and equivalent widths of the PiS line towards 
the observed sources. 



with the large interval in which range the optical depth 
values derived from the fitting procedure. 

The presence of local extinction in the envelopes of the 
dusty sources is consistent with the bow shock model of 
Tanner et al. (2002, 2003). Coteraet al. (1999) had already 
shown that several of these sources are indeed offset from 
nearby local maxima in the extended dust emission and 
temperature distribution. Especially for IRS 21, Tanner 
et al. (2002) indicate that the extended dust emission of 
this source is consistent with a bow shock created by the 
motion of such a massive hot star through the dust and gas 
of the mini-spiral. It is likely that the bow shock scenario 
may be applicable to most of the dust embedded sources 
in the central stellar cluster. For IRS 3, however, the bow- 
shock scenario may not apply. Gezari et al. (1985) find 
IRS 3 as the most compact and (together with IRS 7) 
hottest, bright source (T'^400 K) in the central cluster. 
As suggested in the case of IRS 21 by Tanner et al. (2002) 
IRS 3 may be an optically thick dust shell surrounding a 
mass-losing source, such as a dusty recently formed WC9 
Wolf-Rayet star. 

In addition, a blackbody temperature of ^ lOOOK is 
found for the spectrum representing the highly reddened 
sources located in the North of IRS 13. While it cannot 
be excluded that the individual objects contained in this 
source complex are lower luminosity analogues of the class 
of bow shock objects found by Tanner et al. (2002) and 
Rigaut et al. (2003), their temperature and luminosity is 
well in agreement with the low temperatures of the YSO 
classified by Ishii et al. (1998). 

Detailed modeling, similar to the studies by Tanner 
et al. (2003) and Eckart et al. (2004), based on higher 
angular resolution MIR imaging and spectroscopy (using 
AO or interferometry) is required to unravel the nature of 
the highly extinced sources in the IRS 13N association. 

Finally, a trend of correlation is noticed in the 3.4/iTO 
versus S.Ofim optical depths plot. If a real correlation were 
confirmed between these two features, it would suggest 
that the ISM along the line of sight toward the Galactic 



Center is possibly composed of a mixture of diffuse and 
dense material. Moreover, the plots of the optical depth 
values suggests that part of the 3.4/im feature arises prob- 
ably from the foreground ISM and part of it from the local 
medium associated with the individual sources. 
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